Dynamic receptor team formation can explain the high signal 

transduction gain in E. coli 



Rika Albert * Yu-wen Chiu and Hans G. Othmer 
School of Mathematics, University of Minnesota, Minneapolis, MN 55455 

February 9, 2008 



Abstract 

Evolution has provided many organisms with sophisticated sensory systems that enable them to respond to signals in 
their environment. The response frequently involves alteration in the pattern of movement, either by directed movement, 
a process called taxis, or by altering the speed or frequency of turning, which is called kinesis. Chemokinesis has been 
most thoroughly studied in the peritrichous bacteria Escherichia coli, which has 4 helical flagella distributed over the cell 
surface (Turner et al. 2000), and swims by rotating them (Berg & Anderson 1973). 

When rotated counterclockwise (CCW) the flagella coalesce into a propulsive bundle, producing a relatively straight 
"run", and when rotated clockwise (CW) they fly apart, resulting in a "tumble" which reorients the cell with little translo- 
cation (Berg & Brown 1972, Mcnab & Orston 1977, Turner et al. 2000). A stochastic process generates the runs and 
tumbles, and in a chemoeffector gradient runs that carry the cell in a favorable direction are extended. The cell senses 
spatial gradients as temporal changes in receptor occupancy and changes the probability of CCW rotation (the bias) on 
a fast time scale, but adaptation returns the bias to baseline on a slow time scale, enabling the cell to detect and respond 
to further concentration changes (Block et al. 1982, Segall et al. 1986). The overall structure of the signal transduction 
pathways is well-characterized in E. coli, but important details are still not understood. Only recently has a source of gain 
in the signal transduction network been identified experimentally, and here we present a mathematical model based on 
dynamic assembly of receptor teams that can explain this observation. 

E. coli has five receptor types, but most is known about the aspartate receptor Tar, which communicates with the 
flagellar motors via a phosphorelay sequence involving the CheA, CheY, and CheZ proteins. CheA, a kinase, first au- 
tophosphorylates and then transfers its phosphoryl group to CheY. CCW is the default state in the absence of CheYp, 
which binds to motor proteins and increases CW rotation. Ligand binding to Tar reduces the autophosphorylation rate of 
CheA and the rate of phosphotransfer, and thereby increases the bias. This is the excitation component of the response. 
Bacteria also adapt to constant stimuli, and this is effected by changes in the methylation state of Tar. Tar has four residues 
that are reversibly methylated by a methyltransferase, CheR, and demethylated by a methylesterase, CheB. CheR activity 
is unregulated, whereas CheB, like CheY, is activated by phosphorylation via CheA. Thus, the receptor methylation level 
is regulated by feedback signals from the signaling complex, which can probably shift between two conformational states 
having different rates of CheA autophosphorylation. Attractant binding and demethylation shift the equilibrium toward a 
low activity state of CheA, and attractant release and methylation shift the equilibrium toward a high activity state. These 
key steps, excitation via reduction of the autophosphorylation rate of CheA when Tar is occupied, and adaptation via 
methylation of Tar, have been incorporated in mathematical models of signal transduction (Spiro et al. 1997, Barkai & 
Leibler 1997, Morton-Firth et al. 1999). 

E. coli can sense and adapt to ligand concentrations that range over five orders of magnitude (Bourret et al. 1991). In 
addition, the transduction pathway from an extracellular ligand to the flagellar motor is exquisitely sensitive to chemical 
stimuli. Bacteria can detect a change in occupancy of the aspartate receptor as little as 0.1 — 0.2%, corresponding to 
the binding of one or two ligand molecules per cell. The gain of the system, calculated as the change in rotational bias 
divided by the change in receptor occupancy, was found to be about 55 (Segall et al. 1986), and a long-standing question is 
what the source of this high sensitivity or gain is. Three main sources of gain have been suggested: (i) highly cooperative 
binding of CheYp to the motor proteins, (ii) regulation of CheZ activity , and (iii) indirect activation of many receptors by a 
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ligand-bound receptor. However, it is known that the high sensitivity is present in CheZ mutants (Kim et al. 2001), thereby 
ruling our the second possibility. Furthermore, it was shown that in the absence of cooperativity in signal transduction 
upstream of the motor, a Hill coefficient of at least 11 was needed in the response of the motor to CheYp to explain the 
observed gains of 3 — 6 (Spiro et al. 1997). Cluzel at al. (2000) have confirmed this prediction experimentally, showing 
that the apparent Hill coefficient in the functional dependence of the bias on CheYp is about 10. However, this cannot 
account for all the observed gain, and Sourjik and Berg (2002) have shown, using fluorescence resonance energy transfer, 
that the stage between aspartate binding and CheYp concentration has an amplification 35 times greater than expected. 
None of the existing models of the full signal transduction system (Spiro et al. 1997, Barkai & Leibler 1997, Morton-Firth 
et al. 1999) address this source of gain. 

Receptor interaction, either directly via clustering, or indirectly via an intracellular signal, is a likely source of the 
upstream component of the gain. Receptors are normally dimeric, and it has been suggested that ligand binding affects 
the spatial packing of the receptor array (Levit et al. 1998, Parkinson J. S., University of Utah, personal communication, 
1999). Recent experiments show that transmembrane signaling occurs via receptor clusters or teams, probably of trimers 
of dimers (Ames et al. 2002, Kim et al. 2002). It was previously suggested in analogy with Ising models that clustering 
may enhance the sensitivity at low signals, but it is difficult to obtain both high gain and a wide dynamic range in models 
of this type (Bray et al. 1998). Moreover, while these types of models address the possibility of cooperative interaction as a 
mechanism for generating gain, the nature of this interaction is not specified and thus experimental tests are difficult. More 
recently an abstract model based on the energetics of interactions between receptors was proposed and analyzed by Mello 
& Tu (2003). The model assumes that each receptor dimer can be in an active or inactive state, and that transitions between 
these states are rapid compared to ligand binding. Thus receptors flicker 'on' and 'off between these states, according to 
an equilibrium distribution, and ligand binding biases the proportions in the two states. Parameters can be found so that 
the model reproduces existing data, but again there is no molecular mechanism that can be tested. Our goal here is to 
provide a more mechanistically-based description of the origin of high gain. The model is based on the idea that teams 
of receptor dimers assemble and dis-assemble dynamically, and that different types of receptors can assemble in different 
types of teams. In our analysis assembly and dis-assembly may occur on comparable time scales, but a static scheme in 
which teams exist for long time periods is a limiting case of the model. 

Before describing the model, we observe that the large gain upstream of the motor can be qualitatively understood, once 
the experimentally-determined activity curves are known. The output of the signal transduction network as a function of 
attractant concentration has been studied in several recent experiments, both in vitro (Li & Weis 2000, Bornhorst & Falke 
2001, Levit & Stock 2002) and in vivo (Sourjik & Berg 2002). In vitro experiments use receptor-CheW-CheA complexes 
reconstituted in the presence of attractant and measure the CheA activity immediately after the addition of ATP (Li & 
Weis 2000, Bornhorst & Falke 2001, Levit & Stock 2002). The in vivo experiment of Sourjik & Berg (2002) follows the 
immediate changes in CheYp dephosphorylation after step changes in attractant concentration. These experiments show 
that the measured decrease of the CheA activity with increasing attractant concentration is functionally similar, but not 
identical, to the decrease of the ligand-free receptor concentration. The experimental curves of kinase activity as a function 
of ligand concentration can be fitted with Hill functions of the form 

A(L) = A (l - $) = A (l - k hI l h ) , (1) 

where A represents the measured activity of the network, A$ is the maximal activity in the absence of ligand, ^ is the 
fraction of activity suppressed by ligand binding, L is the ligand concentration, and Ka is the ligand concentration that 
produces half-maximal activity. If we assume that there is no interaction between receptors, the fraction bound with ligand 

is 

e "=K^TL= 1 - e f (2) 

where Kd is the inverse of the affinity for ligand and 0/ is the fraction of receptors free of ligand. If there are only two 
possible states of the receptor complex, free and ligand-bound, and only the former lead to autophosphorylation of CheA 
and a measurable activity, then the activity would have the functional form A(L) = A^6f. However, the experimental ob- 
servations indicate a more complex relationship, in that Ka can be either larger or smaller than Kd and the Hill coefficient 
H can be between 1 and 3 (see Li & Weis 2000, Bornhorst & Falke 2001, Sourjik & Berg 2002, Levit & Stock 2002). 

If the activity is given in the form at Q, we can compute the relative change in activity A(L) and the relative change 
in receptor occupancy for a small change in ligand concentration. Then the gain, which we define as the ratio of relative 
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and this is monotone increasing with L. Thus high amplification is always possible for a sufficiently large ligand concen- 
tration, e.g. L ;» Ka, Kd , and this conclusion holds even if the ligand occupancy has a more complicated dependence 
on L, as long as it approaches one for large L. The explanation of the high amplification is clear from at high ligand 
concentrations the fraction of the activity suppressed tp approaches one, whereas the fraction of receptors free of ligand 
8f approaches zero. Thus the existence of high gain near saturation follows from the functional form of the input-output 
relation of the upstream signal transduction network, and even the simplest assumption of output proportional to 9f leads 
to high amplification for L large compared to Kb- Accordingly, the objective of a model should be to predict the maximal 
activity Aq, the apparent dissociation constant Ka and the Hill coefficient H . 

It is found experimentally that Aq, Ka and H depend on the methylation state of the receptors and the presence 
or absence of the methyltransferase CheR and the methylesterase CheB. Aq increases with the methylation level of the 
receptors and varies approximately 30-fold (Li & Weis 2000, Bornhorst & Falke 2001, Sourjik & Berg 2002). Ka also 
increases with methylation state, and varies over two orders of magnitude (Li & Weis 2000, Bornhorst & Falke 2001, 
Sourjik & Berg 2002, Levit & Stock 2002), which implies that the simplest assumption that ligand-free receptors determine 
the output is not valid (Levit & Stock 2002). The Hill coefficients of the output curves obtained in different experiments 
vary between 1 and 3, and depend very weakly on the methylation level. In vivo experiments also suggest that CheR and 
CheB have a direct effect on the network output, in addition to determining the methylation state of the receptors, because 
the CheR and CheB single mutants show a qualitatively different response than CheRCheB mutants with fixed methylation 
states (Sourjik & Berg 2002). 

Our objective here is to propose a mechanism, based on receptor clustering to form active teams, that can reproduce 
the methylation-induced variability in the network output. There are several recent indications that the receptor-CheW- 
CheA complexes are not static and do not have a one-to-one stoichiometry, as assumed previously. Instead, an oligomer of 
multiple receptor dimers, including different types of receptors, forms the core of an active signaling complex (Ames et al. 
2002, Francis et al. 2002). Since chemotaxis receptors tend to be clustered at one end of a bacterium (Maddock & Shapiro 
1993), we assume that individual homodimers exist in a dynamic equilibrium among singles, teams of two (twofolds), and 
teams of three (threefolds), and that the distribution among st these states depends on the ligand concentration. Our central 
hypothesis is that only threefolds can form complexes with CheW and CheA and activate the autophosphorylation of CheA. 
Because the experimental results we set out to explain all focus on the initial changes in kinase activity, we do not consider 
the slower methyl-transfer reactions. Since phosphotransfer from CheA to CheY is faster than the autophosphorylation of 
CheA, the concentration of phospho-CheY is proportional to the concentration of phospho-CheA, and the output of the 
network is taken to be proportional to the concentration of ligand-free threefolds in the model. 



The basic units of the model are receptor dimers, and we first restrict attention to the inter-dimer association/dissociation 
and the ligand-binding and release reactions for a single receptor type. Homodimers are denoted by R\, twofolds by i?2, 
and threefolds by R3 (see Figure^. Receptor teams can have as many ligand-bound states as there are receptor dimers in 
the team. For example, the R 2 L state contains a single ligand-bound receptor, while R2L 2 has two ligand-bound receptors, 
one on each dimer. We do not consider the state in which two ligand molecules are bound to a homodimer because this 
is energetically unfavorable. We assume that the ligand binding affinity of ligand-free homodimers in a team is the same 
regardless of the binding state of other homodimers in the same team 1 . We allow for the possibility that dimers in receptor 
teams do not have the same affinity for ligand as single receptor dimers {i.e. I2 and ^3 can be different than We assume 
that both ligand-free and ligand-bound receptor dimers can associate to form teams, possibly with different rates i.e. k\ 
and &3 can be either equal or different. 

The main assumption of our model is that ligand binding destabilizes receptor teams and consequently they break into 
smaller units. We allow for every combination of resulting components, but assume that receptors will not release their 
ligands in the process (see Figure 0. In our model only the ligand-free threefolds lead to CheA activation. As a result, 
kinase activity is proportional to the concentration of R3, and its predicted dependence on ligand concentration can be 
compared with the experimental results on kinase activity. 

The kinetic equations for the ligand-free states in Figure [0 are as follows; equations for the remaining states can be 

'The proportionality factors 2 and 3 in the ligand-binding reactions arise from combinatorial effects. 
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Figure 1 : The detailed reaction network for team formation and ligand binding when there is only one type of receptor. 
Individual receptor dimers {R\) can associate to form twofolds (R 2 ) and threefolds (R 3 ). Ligand binding to receptor teams 
leads to the dissociation of the team. Only ligand-free threefolds can initiate kinase activity, 
derived assuming mass-action kinetics. 
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If we define the equilibrium constant K\ 
at L = 0, i?3 satisfies the equation 
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ki/k-i (K 2 = k 2 /k- 2 ), for the formation of twofolds (resp., threefolds); then 
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where Rt is the total receptor concentration fixed at 8fiM. This equation has a unique positive root that tends to zero as 
K\ and/or K 2 tend to zero, approaches its maximum Rt/3 as K\ and/or K 2 tend to infinity, and increases monotonically 
between these limits along any ray in the Ki — K 2 plane. Since the rates of team association/dissociation are not known, 
we assume that K\ — K 2 = K and choose the individual rates k l7 k 2 , fc_ 1; fc_ 2 accordingly. Then R 3 is completely 
determined by K for fixed Rt, and varies with K as shown in Figure|2] 

There is a close parallel between Figure|2]and the experimental observations regarding the dependence of the ligand-free 
output, Ao, on the methylation state of the receptors. The output of the model, R 3 , increases with the team formation con- 
stant K in a very similar way to the increase of Aq with methylation. Notice that the nonlinear increase of R 3 can explain 
both the observation that the lowest methylation state's activity is only a fraction of the activity of the wild type (Sourjik 
& Berg 2002) and the result that the activity of the higher methylation states is very close (Levit & Stock 2002). Based 
on this parallel, we identify different methylation states with different choices for the parameter K, and we choose these 
values such that the ligand-free activities corresponding to these values have approximately the same proportions as the 
experimental measurements of Bornhorst and Falke (Bornhorst & Falke 2001). Thus we identify the unmethylated (EEEE) 
state with K(0) = 10~ 2 \iM~ x and the totally methylated (QQQQ) state with the saturation limit K(A) = 10 3 /iAf~ 1 . 
Bornhorst and Falke constructed 16 engineered states corresponding to all the possible combinations of glutamate (E) and 
glutamine (Q) residues. Based on their results, and similarly to other models (Spiro et al. 1997, Barkai & Leibler 1997, 
Morton-Firth et al. 1999) we assume that the total methylation level is the crucial characteristic of a given state, and not 
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the exact residues that are methylated. We choose the team formation constants for partially methylated states as follows: 
methylation level one (e.g. EEEQ) K(l) = 10~VM -1 , methylation level two (e.g. QEQE) K(2) = 1 ^iM~ l , and 
methylation level three (e.g. QQQE) K(3) = 10 fiM' 1 . 




Figure 2: The dependence of the active team concentration on the team association/dissociation characteristics as cumulated 
into the parameter K. In the following we model different methylation levels by choosing K values such that the ratio of 
their corresponding activities is close to the experimental results of Bornhorst and Falke (2001). 

In order to compare our results with the experimental results in Bornhorst & Falke (2001) and Sourjik & Berg (2002), 
we assume that the receptor is Tar and that the ligand is methyl-aspartate. It is known that the affinity of Tar to methyl- 
aspartate is about ten times less than to aspartate; therefore, we assume that the ligand release rate of a single receptor 
dimer is ten times larger than the release rate of aspartate, which is 70s -1 , while the ligand binding rate is the same as 
the binding rate to aspartate, which is 70/iAf~ 1 s~ 1 (Spiro et al. 1997). Correspondingly, we assume Zi = 70/xM _1 s _1 , 
l-i = 700s -1 . We assume that the ligand release rates of receptor twofolds and threefolds are the same as the release rate 
of a single receptor dimer, i.e. 1-2 = I— 3 = 700s -1 . 

We consider that the association rate of i?i with R2 is slightly smaller than the association between two i?is, and 
the dissociation rate of R3 is proportionally smaller than the dissociation rate of R2, such that the ratios K\ = 
and K~2 = fc2/&-2 are equal. Thus we choose the rates to be fc_i = 0.1s -1 , fc_2 = 0.05s _1 , and we vary k\ and &2 
according to k\ = l0~ 2+d /xM^ 1 s^ 1 , fc 2 = 5 • 10~ 3+d fiM^s^ 1 for d = 0,1,2,4, such that K corresponds to the 
different methylation levels described above. 

As there is no information about the relative rates with which free or ligand bound receptor dimers associate/dissociate, 
we assume that the ligand binding state does not influence team formation and therefore fca = k%, fc_3 = fc_i and 
fcs = &4 = &2, fc_5 = fc_4 = fc_2- For the breakdown of ligand-bound threefolds we consider that processes involving a 
single dissociation are equiprobable, i.e. k$ = kj = kg, = kg = 0.7s -1 , while the processes involving two dissociations 
are less likely, k' 6 = k' r = k' s = 0.07s . 

First we assume that the ligand binding rate of dimers that are part of receptor teams is the same as the ligand binding 
rate of separate receptor dimers, i.e., I2 = h = h- Figure[3]shows the steady-state value of R3 as a function of the ligand 
concentration for four different K values corresponding to four methylation levels. These curves are obtained by solving 
the entire system of steady-state equations using the software package AUTO (Doedel 1981). All curves can be fit with 
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Hill functions of the form ([Q, wherein A(L) = Rs(L). 
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Ligand ( \i M ) 

Figure 3: Concentration of ligand-free threefolds, R3, as a function of external ligand concentration for four differ- 
ent methylation levels, assuming that receptors in teams have the same ligand binding affinity as isolated receptors. 
Circles, K = 1000 /xM" 1 (QQQQ); squares, K = 10 fiM- 1 (QQQE); diamonds, K = (QEQE); triangles, 

K = 0.1 fiM^ 1 (QEEE). The continuous lines represent fits of the form Q- The apparent dissociation constants and Hill 
coefficients are QQQQ - K A = 1.65 fiM, H = 1.37; QQQE - K A = 0.84 fiM, H = 1.23; QEQE - K A = 0.48 fiM, 
H = 1.19; QEEE- K A = 0.30 fiM, H = 1.23. 

This figure captures many features of the experimentally-observed decay in kinase activity for increasing ligand con- 
centrations. The curves are in qualitative agreement with those reported in Li & Weis (2000), Bornhorst & Falke (2001), 
Sourjik & Berg (2002), and Levit& Stock (2002); the apparent K A increases consistently with methylation level while the 
Hill coefficient does not. Thus our scheme provides a possible explanation for the apparent dependence of the receptor 
affinity on the receptor methylation level. We do not change the true affinity for ligand, but varying the affinity of receptor 
dimers for other receptor dimers leads to the differential response of kinase activity to ligand. Unfortunately, the K A val- 
ues predicted by the model are lower than the experimental results obtained for CheRCheB mutants, and the range of their 
variation is also much smaller. Note, however, that the wild type response measured by Sourjik and Berg (2002) exhibits a 
small K A = 1 ijlM, close to our result of K A = 0.48 fiM. 

It is easily seen that the larger K A values observed in experiments on CheRCheB mutants could be explained by 
assuming that receptor teams have a lower affinity for ligand than individual receptor dimers. The cause of this lowered 
affinity could be the close proximity of receptors in teams. To illustrate this case we consider I2 = h = Zi/100, while 
keeping all other parameters at their previous values. The resulting activity curves as a function of the ligand concentration 
are shown in Figure 0] 

These curves agree well with the results given in Figure 2(c) of Bornhorst and Falke (2001). The Hill coefficients 
in (Bornhorst and Falke 2001) range from 1.1 to 2.2, and our values are in this range; the range of K A values are from 
15 fiM (QEEE) to 97 fiM (QQQQ), ours are around this range, too. Moreover, the apparent K A increases dramatically 
with methylation level, while the value of H is about the same for each methylation level; which is consistent with all the 
experimental results (Li & Weis 2000, Bornhorst & Falke 2001, Sourjik & Berg 2002, Levit & Stock 2002). 

Comparing Figure|3]and Figure|4] we can notice that K A depends strongly and inversely on the ligand-binding rates of 
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Figure 4: Concentration of ligand-free threefolds, i?3, as a function of external ligand concentration for four different 
methylation levels, assuming that the affinity of receptors in teams is 1% of that of isolated receptors. The Ka values 
and Hill coefficients are: QQQQ - K A = 112.98 ^M, H = 1.38; QQQE - K A = 37.20 fj,M, H = 1.32; QEQE - 
K A = 16.57 fiM, H = 1.30; QEEE - K A = 6.68 fiM, H = 1.28. 

the receptor teams. Consequently, a possible effect of CheR and/or CheB that would modify the ligand-binding affinity of 
receptor teams would explain the methylation-independent variation of Ka- 

In order to better understand the dynamics of the reaction network, we consider the changes in the concentrations of 
different states on Figure[2under changes in parameters and the ligand concentration. Because the total amount of receptor 
is fixed, changes in ligand propagate through the network until a new steady state is reached. In Table [2 we compare 
the concentrations of the ligand-free and completely occupied states for two external ligand concentrations, two different 
methylation levels, and two types of team ligand binding behaviors. 

As we saw earlier, at L = the proportions of concentrations in these states depends strongly on the methylation state; 
in the highly methylated case the vast majority of receptors are in threefolds, while in the wild type methylation case the 
states are equilibrated. The addition of L = 100 fj,M induces dramatic changes in the state occupancies. These changes 
depend both on the methylation level and the ligand affinity of receptor teams. For h = fa = fa, ligand-free teams all 
but disappear, in both methylation states. The majority of the receptors are now in the ligand-occupied states. In the high 
methylation case the totally ligand-bound teams are most abundant, while in the wild type methylation case almost every 
receptor is in the ligand-bound single dimer state. For fa — fa — h/100, the effect of ligand is much weaker, and it 
depends crucially on the methylation state. At high methylation, the ligand-bound states are sparsely populated, and the 
majority of receptors is in the ligand-free teams, R2 and R3. Notice that the concentration of Tfa almost quadruples, and 
the moderate value of R\L is the only indication of the presence of ligand. In the wild type methylation state the most 
occupied states are the ligand-bound isolated receptors and, still, ligand-free intermediary teams. We can thus conclude 
that in the fa = fa = fa case the most important response to ligand is a vertical flow from ligand-free to ligand-bound 
states, while in the fa = fa = Zi/100 case the most important flow is a horizontal one from receptor teams to individual 
receptor dimers. 

The complexity of the state space induced by team formation also raises the question whether the fraction of ligand- 
bound states is dependent on the methylation level, and how it compares to the ligand affinity of individual receptor dimers. 
We have calculated the dependence of ligand occupancy, (R\L + 2R2L + 3 R3L + 2R2L 2 + 3R^L 2 + 3R 3 L 3 ) / Rt on 
L, and find that it depends on the assumptions about team affinity and methylation state. For example, we find that the 
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apparent dissociation constant corresponding to wild type methylation in the h = h = h case is slightly smaller than 
the Kd — 10 fiM of individual receptor dimers, while in the I2 = I3 = Zi/100 case the apparent dissociation constant is 
higher than Ko- 
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Table 1: The effect of changes in parameters and ligand levels on the distribution of states in the network. K = IfiM 
corresponds to QEQE (wild type), whereas K = 10 3 /iAf corresponds to QQQQ. All concentrations are measured in \iM. 

Finally, we study how the different rates ki and U change the activity curves of the model. 

1. In our choice of parameters we assumed that the ligand binding state does not influence the receptor association 
rate. To explore the effects of different rates, we test the extreme situation when one sets the rates fca and fc_3, 
&4 and fc_4 or fcs and fc_5 to zero with the other parameters and conditions fixed. We find that for each case individual 
concentrations were altered but there was little effect on R3. This suggests that the equality of these rates is not a 
strict condition for the success of the model, and the association between ligand free receptors has the dominant 
effect on R3 

2. In our model we assumed that the ligand-free threefolds constitute the kinase-activating state. To test whether our 
conclusions are generally valid for teams comprised of different numbers of receptor dimers, we assume that the 
association of twofolds with individual dimers is prohibited, and ligand-free twofolds are the kinase-activating state. 
In other words, we set fc 2 = k^ = k 5 = 0. We find that R 2 (L = 0) follows a curve very similar to Figure|2] with the 
only difference that the saturation value for high K is Rq/2 instead of Rq/3. Selecting the same K values for the 
different methylation levels as before we obtain that Ka varies between 2.22/iM (QEEE) and 3.44/iM (QQQQ) 
for l\ = I2 and in the range 25.3/1 M (QEEE) - 314/xM (QQQQ) for I2 = ii/100. The closeness of these results 
to our original model suggests that the number of steps involved in kinase-activating team formation does not have a 
crucial role. 

3. The effect of the single receptor ligand-binding rate l\ on the activity curve is not as strong as the ligand-binding rates 
for receptor teams. When we assume l\ =h = h = 0.7/xM _1 s _1 ratherthan l\ = 70 (j,Al~ 1 s~ 1 , 12 = h = ii/100 
as in Figure 0] Ka increases considerably, but not as much as the change between Figure [5] and Figure |4] This 
suggests that an overall less-than-expected affinity to methyl-aspartate might be at the root of the large observed Ka 
values. 

In conclusion, our results show that a model based on active threefolds of a pure receptor can explain the in vitro 
activity curves (Li & Weis 2000, Bornhorst & Falke 2001, Levit & Stock 2002). The assumptions of methylation depen- 
dent dynamic team formation and ligand-induced breakdown lead to differential kinase activity curves without invoking 
methylation-induced changes in ligand affinity. We turn next to the experimental observations on mixed receptor types 
(Sourjik and Berg 2002). 

Mixed Receptor Types 

Bacteria have several types of receptors, and it is possible that different types of receptor interact in order to be able to 
respond optimally to diverse environmental stimuli. Indeed, the experiments of Sourjik and Berg (2002) suggest that under 
certain conditions both the Tar and Tsr receptors respond to methyl-aspartate. According to these experiments, CheRCheB 
mutants with fixed methylation levels have two apparent dissociation constants corresponding to the Tar and Tsr receptors, 
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respectively, and can be fit by Hill functions of the form 



The simplest model suggested by these results is based on the assumption that the output of the composite system is the 
sum of two individual outputs similar to ([0, one for each of the two pure receptor populations. For this model the total 
output is 

ML) _ A T0 L H ? A SQ L H * 

A A Q L h t + K% T A L Hs +K% S ' 

where A = Ato + Ag . If the Tsr methylation level is constant, as it appears to be experimentally (Sourjik & Berg 
2002), Equation Q predicts that an increase in the Tar methylation level leads to an increase in j3 = Ato/(Ato + Ago). 
However, the experiments indicate that for the CheRCheB mutants, j3 decreases from 0.65 to 0.27 as the methylation state 
of the Tar receptor changes from EEEE to QQEQ. Moreover, (3 appears to be 1 for CheR mutants that are in the lowest 
methylation state and it is for the CheB mutants that are in the highest methylation state (Sourjik & Berg 2002). Therefore 
the experimental results cannot be explained if it is assumed that receptors act independently, which strongly suggests that 
there are interactions between different receptor types, in addition to the interactions within pure types. This leads to the 
possible formation of mixed teams, and in the following we determine whether the model for team formation of pure types, 
extended to two types of receptor, is able to generate response curves similar to those in (Sourjik & Berg 2002). 

We denote the two types of receptors by R and P, and assume that they have different affinities for ligand. We also 
assume that two receptor dimers can associate to form pure or mixed receptor teams, and in the general case the asso- 
ciation/dissociation constants of two R (resp, P or R and P) receptor dimers, K = k\/k-x, (resp., H = h\/h-x and 
M = mx/rri-i), are different. 

The full scheme of all mixed twofold and threefold states and the transitions between them includes 29 states instead 
of the 9 in Figure and involves 38 unknown rates for the receptor association/dissociation reactions alone. It is not 
worthwhile to tackle this level of computational complexity in the absence of any experimental information, and our 
previous analysis suggested that a reduced scheme with receptor twofolds as the ligand-activating state leads to similar 
results as the original scheme. Thus, to reduce the complexity of the analysis of mixed team formation, we only consider 
receptor twofold formation, and we do not consider teams with more than one ligand-bound receptor. Consequently, the 
reaction scheme contains two types of receptor dimers, R\ and Pi, the ligand-bound states of these dimers, three different 
twofolds, i?2, R\Pi and Pi, and six ligand-bound states of these teams, R\L, R2L, R\LP\, R\P\L, P\L and P2L. As in 
the case of pure receptors, we allow for the possibility that teams have lower ligand affinity than homodimers. We assume 
that binding of ligand to either dimer in a twofold induces the dissociation of the twofold, and that the rate of dissociation 
is the same for all teams. Figure[5]shows the entire kinetic scheme and the associated rates. 

This scheme leads to the following steady state equations. 

2fci(l + a 2 L) , 2m 1 \l + (a 2 + b 2 )L] , , , 

K-i + K-2C12L m_i + fc-2(a2 + 02)L 

* hli lt b fr P ? + 2mi[ lt { 7 + !t^ R ^ + C 1 + ^ - pt W 

ft_i+fc_202^ m_i + K_ 2 (a2 + b 2 )L 

where a\ = h/l-i is the affinity of a single R dimer for ligand, 61 = ii/i-i is the affinity of a single P dimer, 
0.2 = h/l-2 is the affinity of an R dimer that is part of a (pure or mixed) twofold, and 62 = « 2 /«- 2 is the affinity of a P 
dimer that is part of a twofold. 

We assume that the kinase-activating output of this system is the concentration R2 + R\Pi + P2 of free twofolds. 
Expressing each of these terms as a function of single receptor concentration we obtain 

m= 1 2m ^ 1 h ^ (10 ) 

fc-i + k-2a%L m-i + fc_ 2 (a2 + b 2 )L h-\ + h- 2 k2L 

The main differences between this output and the simple assumption of non-interacting receptors are the existence of 
the second term depending on R\P\, and the fact that the steady state concentrations of Ri and Pi are coupled. 

Equations (|8} and (|9} can be solved numerically to obtain the output of the network as a function of the ligand con- 
centration. To account for the results reported in (Sourjik & Berg 2002), we assume that the P receptor corresponds to Tsr 
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Figure 5: The kinetic scheme for formation of pure and mixed receptor twofolds. Both receptor types bind ligand, but with 
different affinities, and binding induces team dissociation. 

and its affinity for methyl-aspartate is 10 3 times lower than the affinity of the R receptor (Tar). The other rates are chosen 
to correspond with the rates used in the pure population. Thus we set ai =0.1 /1M" 1 , b\ = 10~ 4 fiM , fc_i = h-i = 
m-i = 0.1 8~\ k- 2 = 70 a -1 , and we allow fci to vary between 10~ 3 i /xM" 1 s" 1 (EEEE) and 10 2 /xJVf _1 s _1 (QQQQ). 

To model the wild type activity curve, we assume that Tar is in its QEQE methylation state, and set fci =0.1 /iM _1 s _1 . 
To capture the surprisingly fast decay of the wild type activity, we assume that receptor teams have the same ligand binding 
affinity as single receptor dimers, i.e., a 2 = ai and 62 = b\. We also assume that the association rate of Tsr into pure 
Tsr teams is lower than the association rate of Tar, i.e., hi = fci/100. These assumptions lead to a Hill function of type 
^with a low Ka, in good agreement with the experimental results Sourjik & Berg 2002 (see filled diamonds in Figure 
[6j. In the CheR mutant both receptors are in their lowest methylation levels since they lack the methylating enzyme but 
have the demetylating enzyme CheB. Again, we assume that the association rate of Tsr is lower than that of Tar, i.e., 
k\ = mi = 10~ 3 fiM^s^ 1 (EEEE) and hi = 10~ 5 /iM _1 s _1 . The decrease in Tar methylation state induces the 
decrease of both the ligand-free activity and the apparent dissociation constant; however, the experiments indicate that the 
Ka of the CheR mutant is close to the Ka of the wild type curve. We are able to reproduce this result by assuming that 
teams have a slightly lower ligand affinity than single receptors, i.e. 02 = ai/10 and 62 = 61/IO (see filled triangles in 
Figure |6). 

Next we consider the CheRCheB mutants, and assume that 02 = a\ /100, 62 = &i/100. We find that for the majority of 
choices for fci, hi and mi the output curves can be fit by generalized Hill functions like Equation[6]with two fast-decaying 
regions characterized by apparent dissociation constants that are several orders of magnitude apart. We identify the lower 
dissociation constant, Kt, with the Tar receptor, and the higher, K $, with the Tsr receptor. 

The experimental results indicate that both the ligand-free output and the apparent dissociation constants Kt and K$ 
increase with increasing Tar methylation levels. Additionally, the parameter (3, indicating the relative weight of the Tar 
receptors in the output, decreases with increasing Tar methylation levels. Our model results in a good agreement with these 
conclusions if we assume that interaction between the two receptor types leads to a moderate variability of the Tsr team 
formation rate. The curves marked by open symbols on Figure [6] present our results for four different sets of pure/mixed 
team formation rates. Different curves have Tar association rates corresponding to different methylation states from EEEE 
to QQQE. We assume that the Tsr-Tsr association rates are lower than the Tar and Tar-Tsr association rates, and they also 
increase with Tar methylation, but with a slower rate 2 . We find that under these assumptions our results are in excellent 
qualitative and good quantitative agreement with the CheRCheB mutant results of Sourjik and Berg (2002). 

In a CheB mutant the receptors are in the highest methylation state, since the action of the methylating enzyme CheR is 
not balanced by CheB. We assume that Tar has a very high association rate, both in pure and mixed teams, while the Tsr-Tsr 
association rate is somewhat lower. As in the case of CheRCheB mutants, we assume that 0,2 — cti/100 and 62 — &i/100. 
We find that the output of such mixture has a single apparent dissociation constant in the milimolar range (see filled circles 
in Figure |6j. 
Discussion 

We have shown that the high upstream sensitivity of the signal transduction network is caused by the negative regulation 
between ligand occupancy of the receptors and kinase activity. Since kinase activity decreases with increasing ligand 
occupancy, at sufficiently high attractant concentrations the relative change in kinase activity is much larger than the 
relative change in occupancy. A related general argument indicates that the sensitivity of the signal transduction network, 

2 Our studies indicate that if we keep the Tsr team association/dissociation rates constant, the high-ligand tails of the activity curves coincide, resulting 
in normalized activity curves that have a reversed order compared to the experimental curves, thus the only way to capture the systematic upwards shift 
with methylation is by assuming a variation of the Tsr methylation levels. 
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ligand concentration ((J-M) 



Figure 6: Activity of modeled Tar-Tsr mixtures as a function of ligand concentration. The Tsr receptors are assumed to be 
twice as abundant as the Tar receptors, and their affinity to methyl-aspartate 10 3 less. Open symbols stand for CheRCheB 
mutants, and have a 2 = a x /100 and b 2 = &i/100. Open circles, K = k\/k-i = 10 ^Al -1 (QQQE), H = h±/h-i = 
6/zM" 1 , M = mi/m_i = 10 pM -1 . Open squares, K = 1 /xM" 1 (QEQE), H = 0.5 pM' 1 , M = l^M~ l . 
Open diamonds, K = 0.1 juM" 1 (QEEE) H = 0.05 fiM' 1 , M = 0.1 plM' 1 . Open triangles, K = 10~ 2 ^Af" 1 (EEEE), 
H = 0.005 fiAl -1 , M = 10~ 2 [iM . The curves can be fit by Hill functions like Equation l|6} with two transition regions. 
QQQE-/3 = 0.36, K T = 22.5 nM,H T = 1.38, K s = 18.2 mM, H = 1.76. QEQE - (3 = 0AI,K T = 9M/j,M,H t = 
1.21, K s = 6.5 mM, H = 1.57. QEEE - (3 = 0.5, K T = 3.68 /xM, ff T = 1.13, ifs = 2.4 mM, i? = 1.43. EEEE -, 
/3 = 0.65, if T = 1.67 ^iM, i? T = 1-09, K s = 1.3 mA/, i? = 1.35. Filled diamonds (wild type), K = M = l^M' 1 , 
H = 10- 4 /iAf- 1 ,a 2 = ffli, 62 = h. Filled triangles (CheR mutant), K = M = lO^/xAf" 1 , iJ = 10- 4 ^Af- 1 , 
a 2 = al/10, 62 = 61/IO. Filled circles (CheB mutant), K = M = 10 6 /iM -1 , i? = 100/iAf-\ a 2 = ai/100and 
62 = 61/IOO. These curves can be fit with single Hill functions like Equation[2 The Ka values and Hill coefficients are 
wild type - K A = 0.43 fxM, H = 0.82; CheR mutant - K A = 0.22 fiM, H = LOT, CheB mutant - Ka = 40 mM, 
# = 1.1. The ratio of amplitudes is wt:CheR:EEEE:QEEE:QEQE:QQQE:CheB=l:0.26:0.35:0.94:1.35:1.52:1.6. 
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defined as the relative change in kinase activity in response to a certain percentage change in ligand concentration, depends 
only on the kinase activity suppressed by ligand 



s =- H ^=- H Krn?> 

and consequently approaches — H at L ^> Ka- This implies that the marked differences between the sensitivity of the 
wild type and CheRCheB mutants found by Sourjik and Berg (2002) are caused by the fact that for the ambient ligand 
concentrations studied the wild type receptors have reached the maximum sensitivity, while the CheRCheB mutants have 
not. 

To illustrate this point we include a figure depicting the absolute value of the sensitivity of two receptor populations 
to a 10% change in ligand concentration (similar to Figure 3b in Sourjik and Berg 2002). The first population's kinase 
response is described by a Hill function like Equation^with the parameters Ka — 1/iM, H = 1 while the second's kinase 
response follows Equation[6]with Kt = 150/xAf , Ks = 100mA/ and Ht = Hs = 1- In the ambient ligand concentration 
range 10^M < L < 10 4 ^M the sensitivity of the first population is constant since L > Ka, while the second population, 
having L < Ks, has a varying and much smaller sensitivity. This behavior is in excellent qualitative agreement with the 
experimental observations of Sourjik and Berg (2002). 



10" 



a 



ti -l 
glO 

ID 
00 



[] 



10"? V 



10 



o o o ooooooooQoo6ocsmDoooo6ooQosooo 

K A < range of L 



OD 



□ 



□ 



K > range of L 



10 10 10 

Ligand concentration (|0,M) 



10 



Figure 7: Sensitivity of two receptor populations to a 10% increase in ligand concentration. The sensitivity is defined as 
the ratio of the relative change in kinase activity and the relative change in ligand concentration. 

We have demonstrated that a model based on the assumption that homodimers of a receptor must aggregate into teams 
of three in order to activate the autophosphorylation of CheA can adequately explain the observed dependence of the 
kinase activity on the ligand concentration for a pure receptor. Our model is in qualitative agreement with the experimental 
results, and shows that methylation-dependent kinase activity does not necessarily imply methylation-dependent ligand 
affinity. We also showed that the concentration corresponding to half-maximal kinase activity need not coincide with the 
apparent ligand dissociation constant of the receptor population, nor does the latter coincide with the dissociation constant 
of an isolated homodimer. 

Our model assumes that receptor populations possess a dynamic balance between homodimer, twofold and threefold 
states, as opposed to an ordered threefold structure. This prediction, along with our assumptions for team formation and 
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dissociation rates could be tested experimentally in in vitro receptor preparations. Furthermore, in order to quantitatively 
reproduce the experimental results on CheRCheB mutants (Bornhorst & Falke 2001) within the framework of the detailed 
model it is necessary that twofolds and threefolds of receptors have a lower affinity for ligand than an isolated homodimer. 
This theoretical prediction could be verified experimentally by testing the affinity of homogeneous receptor preparations 
(i.e. only dimers or only teams). 

When there are multiple receptor types, the experimentally-determined activity curves display complex dependence 
on the ligand concentration, but they can be satisfactorily reproduced by our model. One consistent assumption that was 
needed is that the association rate of Tsr teams is lower than the association rate of Tar and Tar-Tsr teams. This assumption 
was vital in reproducing the wild type, CheR and CheB mutant curves, and it suggests the existence of receptor specificity 
in team formation capabilities. This feature could be caused by receptor-specific methyl-accepting activities that were 
confirmed experimentally (Barnakov et al. 1998). 

Our results also confirm earlier suggestions that changes in methylation state (or association/dissociation rates) alone 
cannot explain the qualitative difference between the wild type and CheRCheB mutant activity curves. We were able to 
reproduce the shift by assuming that in CheR or CheB mutants receptor teams have lower affinity for ligand than individual 
receptor dimers. Note that this effect is weaker in CheR mutants, but still existent. 

Our analysis deals only with the early response to changes in ligand concentration, since we have neglected methylation 
of receptors and downstream phosphotransfer reactions. It remains to integrate the model for the early response developed 
here with a complete model such as given in (Spiro et al. 1997) for later events. It is of course feasible to do this 
computationally, but given the complexity of the association scheme for the formation of signaling teams shown in Figure 
[2 it may be difficult to extract qualitative insights analytically. Some simplification exploiting the disparity in time scales 
of the various processes will certainly be needed. 

This work was supported by NIH Grant #GM-29123 to H. G. Othmer. We thank Sandy Parkinson for helpful discussions 
at various stages of the model development. 
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